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Abstract 

We present evidence for a phase transition in a theory of 2D causal set 
quantum gravity which contains a dimensionless non-locahty parameter 
e £ (0, 1]. The transition is between a continuum phase and a crystaUine 
phase, characterised by a set of covariant observables. For a fixed size of 
the causal set the transition temperature /3~^ decreases monotonotically 
with e. The line of phase transitions in the v/s e plane asymptotes 
to the infinite temperature axis, suggesting that the continuum phase 
survives the analytic continuation. 

Causal set theory(CST) is a discrete approach to quantum gravity which 
combines local Lorentz invariance with a fundamental discreteness [T]. The 
spacctime continuum is replaced by a locally finite poset or causal set, with the 
order relation -< being the analog of the spacetime causal order. The continuum 
quantum gravity path integral is thus replaced by a discrete sum over causal 
sets 

zcsT = E cxp^^[^l/^ (1) 

ceo 

where f2 is a sample space of causal sets and S[C] is an appropriately chosen 
action. The combination of local Lorentz invariance with fundamental discrete- 
ness gives rise to a non-locality in the continuum approximation, making the 
extraction of local geometric data highly non-trivial. The recent construction of 
a discrete Einstein-Hilbert action, the Benincasa-Dowker action for causal sets 
[2], thus allows us for the first time to begin a serious study of the causal set 
partition function Zcst- 

Apart from a choice of action, Zcst also depends crucially on the sam- 
ple space ft. A natural starting choice for O is the collection of countable 
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causal sets; in classical sequential growth models of causal sets, for example, ft 
is further restricted to causal sets that are past finite [31 H] . The collection of 
TV-element causal sets rijv is known to be strongly dominated by the "Klcitman- 
Rothschild" (KR) class of causal sets in the large N limit. These arc of a fixed 
time extent with only three "moments of time" , and admit no continuum ap- 
proximation [5]. This presents a potential "entropy problem" in CST. In order 
to be able to recover spacetime-like behaviour, therefore, the causal set action, 
or more generally the choice of dynamics, should be able to counter this en- 
tropy. Indeed, classical sequential growth dynamics is an example in which the 
KR entropy is made sub-dominant by the dynamics [3]. 

In this work we consider a two dimensional theory of causal sets, defined 
by the 2-dimcnsional Benincasa-Dowker action S2d 12, and an order theoretic 
dimensional restriction of fijv to ft2D , the set of A^-element "2D orders" , defined 
as follows. Let S — (1, . . . , N) and U = (wi, ii2, • ■ • u^), V ~ (ui, U2, . . . ww), 
with Ui, Vi € S, Ui ^ Uj, Vi ^ Vj for i ^ j. U and V arc then total orders with 
-< given by the natural ordering < in S: for every pair i ^ j either Ui < Uj 
or Uj < Ui, and similarly for V. An A^-element 2D order is the intersection 
C ^ U nV oi two total A^-element orders U and V, i.e., e,; -< ej in C iff Ui < Uj 
and Vi < Vj. A useful example of an iV-element 2D order is a set of N events 
{ei, . . . , Bat} in 2d Minkowski spacetime ordered by causality, and such that in 
light cone coordinates e; = (ui,Vi), Ui ^ Uj and Vi ^ Vj for i ^ j- 

The motivation for restricting the sample space to fl2D stems from the fact 
that the causal set discretisation of a conformally fiat, topologically trivial 2d 
spacetime is a 2D order [6]. However, not all 2D orders can be approximated 
by continuum spacetimes. This means that the choice of fl2D corresponds only 
to a restriction of poset dimension and not spacetime dimension. Moreover, 
only in a very limited sense do these dimensions coincide: every 2D order C 
admits an order preserving embedding $ into a patch of 2d Minkowski spacetime 
(^M, ry), i.e., for every -< Cj in C, $(6-;) causally precedes $(ej) in i^M^-q). 
Such an embedding though necessary, is not sufficient to ensure a continuum 
approximation for C. 

A striking feature of 2D orders is that in the asymptotic limit N ^ oo, il2D 
is dominated by "random" 2D orders, namely those approximated by Minkowski 
spacetime [6l[7l|8]. A random 2D order is the intersection of two total orders U = 
. . . , un) and V ~ (wi, . . . vn) which are chosen randomly and independently 
from S. Hence unlike the unrestricted sample space, spacetime like causal sets 
dominate the uniform measure on fl2D- It is thus of obvious interest to study 
the effect of 824 on this entropic feature of fl2D- 
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While there is no natural Planck scale in 2D gravity, in CST one requires 
a volume cut-off Vp = Ip in order to realise the continuum approximation, and 
this plays the role of a fundamental scale. In addition, the 2d Benincasa-Dowker 
action S2d[C] for a causal set C includes a "non- locality" scale Ik > Ip required 
to suppress large fluctuations about the mean in the continuum approximation 
[11 E] . S2d [C] can be expressed in terms of the abundances Nn of the "intervals" 
= {k\i -< k -< j} of fixed cardinality 71 in C 

-52d[C,e]=4e(7V-2e^iV„/(n,e)), (2) 

where ( = 1^/11 and 

Figure [1] shows the typical behaviour of /(n, e). 




Figure 1: The function f{n, e) for e — 0.12. 

For n > 1/e, f{n,e) ^ 0, and thus, effectively, the action for a fixed e 
includes the abundance of intervals only upto size 1/e. To avoid infrared errors 
in counting such intervals for finite N, e should in addition be bounded below 
by Ne > 1. 

The goal of writing down a partition function for quantum gravity is to even- 
tually construct and calculate expectation values of covariant observables. The 
analog of covariance in causal set theory is label independence, and hence we 
seek to construct label invariant observables. We have already encountered an 
example of such an observable, namely, the action S2d/f^- Another important 
covariant observable is the Myrhcim-Myer dimension [lOj for a causal set, which 
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in the continuum approximation reproduces the Minkowski spacetime dimen- 
sion. We will construct several such covariant observables in what follows and 
use them to determine the existence of a continuum approximation. 

However, as we have defined it, the sample space is itself the collection 
of labelled 2D orders. While every 2D order admits A^! worth of relabellings, 
each of these relabellings does not produce a distinct labelled 2D order. If C is 
a labelled causal set containing a pair of elements and Cj with the same past 
and future sets, the relabelling z f-^ j is an automorphism Hi^^j : C ^ C . An 
extreme example is the antichain or completely unordered set - every relabelling 
produces the same labelled causal set. Since the action itself is label invariant, 
the partition function over unlabelled causal sets has an additional weight of 
N\/\Aut{C) \ for each unlabelled 2D order C, where |^ui(C)| is the cardinality of 
the Automorphism group of a labelled counterpart. While it is possible to view 
this additional measure as a "quantisation ambiguity", it should be stressed 
that the observables that wc construct arc nevertheless strictly covariant. 

An approach to evaluating the expectation values of covariant observables 
is to Euclideanise the partition function, in the process converting the quan- 
tum system into a thermodynamic one. Replacing the set of Lorentzian metrics 
wholesale for Euclidean ones, however, makes little of the importance of causal 
structure, and moreover, can lead to highly fractal, non-manifold like behaviour 
[llj . Instead, as suggested in jl2j one can transition to a thermodynamic par- 
tition function by introducing a new parameter /3 into the action and taking 
P — > i/3. This leaves the sample space unchanged and hence there is no ambi- 
guity in the interpretation of the covariant thermodynamic results. A similar 
"parameter-based" analytic continuation employed in the causal dynamical tri- 
angulation approach, on the other hand, does lead to a Euclideanisation of the 
sample space itself [13]. 

For concreteness, we adopt the following prescription. Namely, we replace 
the complex weights exp''^^''/'' with exp''^'^^''/'', and obtain the thermodynamic 
partition function by taking /3 — >■ 

Zn= exp-'^^-I^l/^ (4) 

The /3 — >■ limit is the uniform distribution, dominated by the 2D random 
orders, approximated by an interval in Minkowski spacetime for finite N . As 
[3 ^ oo, since the action is not positive definite, causal sets with the largest 
negative values of the action should dominate, modulo cntropic effects. Thus, 
one expects a cross-over at finite /3. 

In this work, we use Markov Chain Monte Carlo(MCMC) methods to study 
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this partition function. The results we present are the first in a larger effort 
to study causal set quantum dynamics using MCMC techniques [ij. The re- 
striction to 2 dimensions leads to a substantial simplification which translates 
into rapid mixing or thermalisation of the Markov Chain. Our simulations are 
carried out for relatively small causal sets {N = 50) but this is sufficient to show 
emergent continuum behaviour. One of the main observations of our present 
work is the existence of a phase transition at finite rather than the cross-over 
suggested above. The transition is from a continuum-like phase to a crysalline 
non-continuum phase, and is well characterised by the change in the expectation 
values of a set of covariant observables. 

We first briefly review the basics of causal set theory, and refer the reader to 
[T3] for more detail. A causal set is a locally finite partially ordered set (posct), 
i.e., a countable collection of elements with an order relation -< which is (i) 
transitive {x -< y ,y < z ^ x ^ z), (ii) irreHexive, {x yi. x) and (iii) locally 
finite, i.e., if Past{x) = {w E C\w -< x} and Fut{x) = G C|z/; >- x} then 
the cardinality of the set Past{x) n Fut{y) is finite. We say that two elements 
are linked if x ^ y and there is no z such that x < z ^ y. Local finiteness 
implements the physical requirement of a fundamental spacetime discreteness: 
a finite spacetime volume contains only a finite number of "spacetime atoms" . 
A continuum spacetime {M, g) is said to be an approximation to an underlying 
discrete causal set C for a spacetime volume "cut-off" Vc, if there exists an 
embedding $ : C — (Af , g) which is (i) order preserving: for x,y E C, x -< y 4^ 
^{x) -<M ^{y), and (ii) $(C) is a Poisson distribution of events in {M,g): 

Pv{n) = ^^^exp-^ (5) 

is the probability of finding n elements of $((7) in a spacetime volume V. Con- 
versely, one can generate a causal set C from {M, g) via a Poisson sprinkling, 
with -< induced by ~<m for the sprinkled points. A causal set generated this 
way is therefore not a regular lattice but a "random lattice" . Starting with a 
sample space of finite element causal sets ^In, the thermodynamic limit N ^ oo 
does not correspond to the continuum limit of the theory, but rather its infrared 
limit. This means that a finite cardinality causal set can be well approximated 
by a region of continuum spacetime without taking the TV — >■ cxj limit. 

We construct the MCMC via the exchange move on a 2D order C = U DV 
defined as follows. First pick U or V at random. Wlog, let this be U. Next, pick 
a pair (ui,Uj) of elements of U at random and perform the exchange Ui <-> Uj, 
while leaving V unchanged. The new 2D order then has the new elements 

= = (uj^Vi) and e^- = (u'^v'^) — (ui,Vj), while all other elements 
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remain the same. It is clear that the move is reversible. 

Claim 1. The exchange move on labelled 2D orders has no fixed points in the 
space of labelled 2D orders. 

Proof: Wlog, let u,; o Uj for some i ^ j- There are three cases to examine: (a) 
If Ci -< Cj , then Ui < Uj and Vi < Vj . Under the exchange u'- > u'^ while v'^ < v'^ 
which means that and e'j are now unrelated, (b) Similarly if Ci >- ej, then 
after the exchange the two elements are unrelated, (c) If they are unrelated to 
start with then either (i) Ui < Uj and Vi > Vj, in which case after the exchange 
e'j < e[ or (ii) u,; > Uj and Vi < vj, in which case after the exchange e'j > e[. 
In all three cases, therefore the exhange move makes a non-trivial modification 
to the relationship between the ith and jth elements, thus giving rise to a new, 
distinct labelled causal set. □ 

We employ a Metropolis-Hastings algorithm, accepting a move if the differ- 
ence in the action AS'2d is negative and rejecting a move only if Qx^-P^S2d/h <^ j,^ 
where r is a random number in [0, 1). Since each move depends on a pair of 
elements, we define a sweep to be N{N — 1) /2 moves, with the following observ- 
ables recorded every sweep, (i) The ordering fraction x — 'ir /N{N — 1) where 
r is the number of relations in the causal set with N{N — l)/2 the maximum 
number of relations possible; x is therefore analogous to the filling fraction in 
an Ising model. In 2 spacetime dimensions x is also the inverse of the Myrheim- 
Myer dimension. Thus, for flat spacetime, we should expect (x) ~ 1/2. (ii) The 
action S2d[C]/h: for flat space it was shown in [16] that {824 1 fi) ~ 4 for e = 1 
and simulations show that this is also true for e 7^ 1. (iii) Time Asymmetry: 
given that the action is itself time-reversal invariant, it is useful to check whether 
the dynamics preserves this property. A rough measure of time asymmetry tas 
is the difference in the number of maximal and minimal elements in the causal 
set; in flat spacetime {tas) ~ 0. (iv) The height h of the causal set or the 
length of the longest chain. This corresponds in the continuum approximation 
to the maximum proper time in flat spacetime [17] . (v) The abundances Nn of 
n-element intervals for all 71 g [0, — 2]. In flat spacetime Nn has a specific 
monotonic fall off with n and this provides a useful comparison. 

Our focus in this work is to examine the behaviour of these observables as 
function of the inverse temperature /3 as one varies the non-locality parameter 
e. We present results for N = 50, with a range of values for e between 0.1 and 1. 
Wc perform the simulations for 10, 000 sweeps which translates to 12.5 million 
attempted moves. 

Wc test for thcrmalisation starting from several different types of initial 2D 
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orders. Our current algorithm includes the following 8 initial types of causal 
sets: (i) a randomly labelled chain or totally ordered set (ii) the antichain (iii) a 
random 2D order, which is approximated by the Minkowski interval (iv-v) two 
2D random orders of size 7V/3 each with an intervening element chain or 
antichain, (vi) two 2D random orders of size N/16 with an intervening antichain 
of size 7N/8 (vii) two iV/2 element random orders stacked on top of each other 
and (viii) a crystalline causal set of the kind we will describe shortly. We find 
that there is rapid mixing starting from these varied configurations, an example 
of which is shown in Figure [5J The expectation values ( O } are calculated from 
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sweeps 

Figure 2: The exchange move gives rise to a rapid thermalisation starting from the 
eight different different starting 2D orders. Here we show a plot of action v/s number 
of sweeps for e = 0.12 and /3 = 0.1. Thermalisation takes place by the 2nd sweep. 

a sampling per autocorrelation time to- As expected, to increases with /3 so 
that the number of independent samples rapidly decreases for larger /3. 

For fixed e we find a rapid change in the expectation values of the obscrvablcs 
around a fixed /3 = /3c suggesting a phase transition as shown in Figure ([3]). 
The exact value of /3c can be determined by looking at the behaviour of the 
autocorrelation times which peak at the same /3c for all observables. For the 
example shown in Figure [3] it is clear that roughly, 2.2 < /3c < 2.5. We have 
included the error bars in the plots, but these are typically very small. 

That there is an actual differentiation into two phases becomes explicit on ex- 
amining the 2D orders themselves. We record the configurations in the Markov 
Chain every 100 sweeps. For /3 < /3c, we find evidence for a "continuum" phase 
(Phase I), and for /3 > Pc, "crystalline" phase (Phase II). We show examples of 
configurations in these two phases in Figure IH where the light cone coordinates 
have been turned clockwise by 7r/4 for ease of plotting. 

A quick look at a typical causal set Cj from Phase I shows a causal set 
that resembles a random discretisation of a spacctime. Indeed, the expectation 
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Figure 3: Observables v/s /3 for iV = 50, e = 0.12. 



values for a fixed j3 < j3c corroborate this. For example for e = 0.12 and j5 ~ 0.1, 
(a) ordering fraction: < x >= 0.499 with an error less than 10"'^ , which means 
that the Myrheim-Myer dimension < d^/A/ 2.004, (b) height: < h >= 
10.232±0.014, which should be compared to the height of a = 50lp Minkowski 
interval which is 10/p, (c) time asymmetry: < tas >= 0.027 ±0.024 and the (d) 
Action: < S > /h ~ 3.846 ± 0.013 which should be compared to the residual or 
boundary value of 4 for Minkowski spacetime [TB]. In addition, if we plot the 
abundance Nn of the intervals of cardinality n in C/ and contrast it that for a 
random 2D order, we find that these match very closely as shown in Figure [SJ 
As /3 nears the transition, the expectation values of these observables gradually 
change. However, the typical causal set continues to retain the features of a 
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Figure 4: On the left is a causal set in the Continuum Phase and on the right one 
in the Crystalline Phase. {U, V) are light cone coordinates that have been rotated 
clockwise by 7r/4. 



random lattice. Thus, looking at the explicit values of the observables, it is 
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Figure 5: The distribution of Nn in C/. This is compared with the distribution for a 
2D Random order which is a discretisation of Minkowski spacetime. 

clear that Phase I corresponds to a continuum phase. 

A typical causal set C// in Phase II on the other hand, has a most unexpected 
character as shown in FigureH) It shares some superficial features of a KR causal 
set, being of limited time extent, but does not strictly belong to this class. It has 
a regularity, or crystalline nature, suggesting that it does not have a continuum 
approximation. The values of most observables differ considerably from those in 
Phase I. For example, for e = 0.12, /3 = 3.5 (a) Ordering Fraction: < x >= 0.579 
which gives a fractal dimension for the causal set of < (Imm >^ 1.727, (b) 
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Height: < h >= 4.180±0.018 (c) Time Asymmetry: < tas >= -10.292±0.102, 
which is a large deviation from the expected value of zero. This may be the 
result of a spontaneous breaking of the time symmetry, similar to that in the 
low temperature phase of the Ising model, but may also be the result of poorer 
statistics at larger /3. (d) Action: < S > /h~ —41. 367 ± 0.054, i.e., the action 
tries to take on the lowest possible (negative) value. In contrast to Phase I, 
the abundances of intervals as shown in Figure [S] is very different from that of 
Minkowski spacetime. In particular, Nn align themselves with the positive part 
of f{n,e) and vanish for those n for which f(n,e) < 0, thus minimising the 
action. 
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Figure 6: The distribution of N„ in Cii. On the left it is compared to a 2D Random 
order for n < 20. On the right it is shown for all n and can be seen to roughly follow 
the positive part of f{n, e). 



The existence of these distinct phases thus strongly suggest a phase transi- 
tion. It is tempting to draw the obvious analogy with the Ising model: at high 
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temperatures one has the disordered or random Phase I, while at low tempera- 
tures, there is the higly ordered crystalhne Phase II which exhibits a spontaneous 
breaking of symmetry. Ahhough it is difheuh at this stage of our work to assess 
the order of this transition, there are hints that it may be of second order. To 
begin with, the autocorrelation time peaks at the phase transition as do the 
fluctuations in the observables. The size of the fluctuations gives us a rough 
estimate of the temperature of the phase transition, and this coincides with the 
peak in the autocorrelation time. Moreover, the transitions shown in Figure [3] 
appear to be smooth. A more conclusive assessment would require a detailed 
understanding of the dependence of these results on the cardinality N. 

It is important to stress that while the nature of this phase transition is 
deflnitely of interest, it does not play a crucial role in determining continuum 
behaviour. In other lattice-based approaches in which discretisation is used 
only a calculational tool, the appearance of a second order phase transition 
signals the fact that the continuum limit of the theory exists. However, in a 
fundamentally discrete theory like causal set theory, it is only the continuum 
approximation we seek, and as we have just seen, this does not depend on the 
existence and the nature of a phase transition. 

More important to our discussion, then is the question of what this ther- 
modynamic calculation can mean for quantum gravity. How much of the above 
discussion, if any, survives the analytic continuation? As e varies in (0, 1], our 
simulations show that the phase transition survives, but the critical tempera- 
ture increases monotonically with e. Using the maximal size of the fluctuations 
to estimate the critical temperature, we plot as a function of e. The negative 

axis corresponds to the region of interest, i.e., the quantum regime, while 
the positive axis corresponds to the thermodynamic regime to which our 
above calculations belong. As shown in Figure [7] the line of phase transitions 
which separates Phase I from Phase II asymptotes to the = axis which 
itself belongs to the continuum Phase I. This strongly suggests that the non- 
continuum crystalline Phase II is confined to the > region and that the 
continuum phase survives the analytic continuation into the < region. The 
phase diagram also suggests that smaller e is in a sense "favoured" by Phase 
I, which is consistent with having smaller fluctuations in the action and hence 
a more reliable continuum approximation. A similar type of analysis has been 
used recently in studying the phase structure of QCD ^TE\ . 

A more careful study of the phase diagram is clearly in order. For example, 
does the gradual change in the observables in Phase I as one moves towards /3c 
signal a significant shift away from Minkowski spacetime? In a recent paper [19j 
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Figure 7: This is a plot of the locus of critical points as a function of e. For 
P'^ > /3c the system is in the crystalline phase and for < /J^ it is in the continuum 
phase. 



the phase structure of the 4d causal dynamical triangulation model of quantum 
gravity was studied. The phase diagram of this theory includes three phases, 
two of which, Phase C and Phase B, seem to be at least superficially analogous 
to our Phases I and II respectively. The phase transition B-C is argued in [T9] 
to be second order and it would be interesting to explore whether these two 
theories lie in the same universality class. 

Since a more ambitious goal is to work with the unrestricted sample space 
n, and an action with at least dimension 4 [14], it is not immediate that 2d 
CST can teach us straightforward lessons. Nevertheless, our analysis opens a 
new window into causal sets, and with it a host of questions that can finally 
begin to be addressed. One of the more interesting of these is whether there is 
a renormalisation group type analysis with stable fixed points for e < 1. This 
would suggest that the non-locality scale is not a free parameter, but can be 
determined from the quantum dynamics. 
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